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Anthony  J.  Midey,  Jr. 

NCEE  Prime  Contract  -0029  from  23  June  1999  to  22  June  2000 

During  the  aforementioned  contract  period,  I  have  employed  a  variable 
temperature-selected  ion  flow  drift  tube  (VT-SIFDT)  and  a  high  temperature  flowing 
afterglow  (HTFA)  to  study  ion-molecule  reactions  from  220-1400  K.  I  have  investigated 
the  effects  of  electronic,  translational,  rotational  and  vibrational  energy  on  the  kinetics  of 
these  reactions.  This  work  has  been  conducted  under  the  direction  of  Dr.  Albert  A. 
Viggiano  at  the  Air  Force  Research  Laboratory  (AFRL)  in  the  Space  Vehicles 
Directorate  (VSBP)  at  Flanscom  AFB,  MA. 

Plasma  chemistry  occurring  at  very  high  temperatures  critically  impacts  Air  Force 
systems.  Plasmas  can  develop  around  re-entry  vehicles  that  interfere  with 
communications.  In  addition,  the  chemistry  of  the  ionosphere  occurs  at  temperatures 
upwards  of  2000  K.  Consequently,  accurate  modeling  of  ionospheric  processes 
requires  measurements  of  reactions  at  high  temperatures. 1  The  HTFA  has  been 
designed^  and  used  to  address  these  issues.'' -  3, 4  However,  weakly  ionized  systems 
are  also  of  interest  and  are  readily  studied  using  the  VT-SIFDT. 

Understanding  high  temperature  ion-molecule  chemistry  also  facilitates  the 
design  of  airbreathing  hypersonic  vehicle  combustors  where  the  ion  chemistry  may 
enhance  hydrocarbon  fuel  combustion. ^  Unfortunately,  ion  chemistry  may  also 
contribute  to  soot  formation  in  hydrocarbon  combustion.®*  ^  Nevertheless,  the  first  step 
in  establishing  this  technology  is  modeling  the  chemistry  occurring.  However,  only 
neutral  molecule  chemistry  has  been  included  in  recent  models  of  soot  generation.® 
Accurate  models  of  these  processes  thus  require  high  temperature  data  for  ion- 
molecule  reactions  with  many  different  hydrocarbons,  including  alkyibenzenes  which 
comprise  about  25%  by  volume  of  standard  jet  fuels.  The  previously  mentioned  HTFA 
apparatus  has  subsequently  been  employed  to  address  these  problems.  Recent 
modifications  to  the  HTFA  have  expanded  its  capabilities,  allowing  measurements  of 
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ion-molecule  branching  ratios  at  temperatures  over  700  K.-*''’*  The  VT-SIFDT  and  the 
HTFA  combined  permit  branching  ratios  and  rate  constants  to  be  measured  from  85  to 
1800  K  at  pressures  from  up  to  2  torr  and  kinetic  energies  up  to  1  eV. 

The  VT-SIFDT  operates  in  the  following  manner. '’2  ions  are  generated  in  a  high 
pressure  source  region  using  electron  impact  on  a  source  gas.  The  desired  reactant  ion 
is  mass  selected  by  a  quadrupole  mass  spectrometer  and  injected  into  the  flow  tube. 
The  ions  are  entrained  in  a  fast  flow  of  helium  buffer  gas  and  thermally  equilibrate  to  the 
flow  tube  temperature,  which  is  controlled  by  a  combination  of  heating  elements  and 
cooling  lines.  After  thermal  equilibration,  the  neutral  reactant  is  introduced.  The 
residual  reactant  ions  and  all  product  ions  are  sampled  with  an  orifice  and  are  mass 
selected  by  another  quadrupole  mass  spectrometer  then  detected.  For  a  given 
temperature,  the  rate  constants  are  obtained  using  the  decline  of  the  reactant  ion  signal 
over  a  previously  measured  reaction  time  as  a  function  of  excess  neutral  reagent 
concentration  where  the  average  center-of-mass  kinetic  energy  is  given  by:''^ 


+  -/c„T 


(1) 


In  Eq.  1 ,  mu  mu  and  mn  are  the  reactant  ion,  buffer  gas  and  reactant  neutral  masses, 
respectively.  The  drift  tube  consists  of  a  series  of  concentric  rings  to  which  a  voltage  is 
applied  to  accelerate  the  reactant  ions  to  a  kinetic  energy  above  the  thermal  value 
described  by  the  second  term  in  Eq.  1 .  The  first  term  above  gives  the  drift  tube 
contribution  to  the  total  kinetic  energy.  Operating  the  drift  tube  with  no  voltage  applied 
allows  the  instrument  to  be  run  as  a  VT-SIFT.  Pseudo-first  order  kinetics  hold  under  all 
of  these  conditions  and  the  rate  constants  are  easily  calculated.  The  branching  ratios 
are  measured  in  two  steps  to  minimize  mass  discrimination  in  the  mass  spectrometer 
and  to  allow  accurate  determination  of  the  individual  reaction  channels  that  can  differ  by 
only  one  mass  unit.  Relative  errors  in  the  rate  constants  are  ±1 5%  and  the  absolute 
errors  are  ±25%,  while  the  relative  errors  in  the  branching  ratios  are  ±25%. '^2 

Rate  constants  and  branching  ratios  are  measured  in  the  HTFA  in  a  similar 


manner,2  with  a  few  exceptions.  The  HTFA  has  no  source  selection.  Consequently, 
branching  ratios  can  be  measured  only  for  ions  that  can  be  generated  with  large  flows 
of  source  gas  because  the  ions  arise  from  interaction  with  high-energy  helium  species 
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formed  in  the  source  region.  These  large  flows  guarantee  that  all  unwanted  reactive 
species  are  consumed  before  the  neutral  reagent  gas  is  introduced.  Corrections  must 
also  be  made  to  the  branching  ratios  for  reactive  species  formed  in  the  flow  tube  of  the 
HTFA  (e.g.,  N/  when  making  N2'^  ions).  The  uncertainties  in  the  HTFA  data^.  9  are 
similar  to  those  in  the  VT-SIFDT. 

The  rate  constants  and  branching  ratios  for  the  reaction  of  H30'^(H20)n  clusters 
with  aldehydes  CH2O  and  CH3CHO  have  been  studied  from  150-294  K  with  n=0-5  and 
200-294  K  with  n=0-4  in  the  VT-SIFDT.  These  measurements  address  assumptions  by 
Yu  et  al.'^4  that  all  of  the  aldehydes  that  have  been  used  to  model  ultrafine  aerosol 
formation  in  jet  exhaust.from  the  interaction  of  charged  particles  with  volatile  organic 
compounds.  The  current  experiments  also  improve  upon  previous  flowing  afterglow 
measurements  where  both  H2O  and  the  aldehyde  are  in  the  flow  tube  which  means  that 
the  system  is  in  equilibrium.'’ 9.  is 

The  current  results  are  given  in  Tables  1  and  2  for  CH3CHO  and  CI-I2O.  The 
CH3CHO  reactions  proceed  at  the  collision  rate  at  both  temperatures  for  n=0-4  with 
ligand  switching  predominating  for  n=1-4.  The  bare  H30‘^  ion  reacts  via  proton  transfer. 
Higher  order  ligand-switching  products  from  secondary  reactions  are  also  observed. 
The  CH2O  reactions  show  similar  products,  but  the  rate  constants  are  both  size  and 
temperature  dependent.  H30^  again  reacts  by  proton  transfer  at  the  collision  rate  at  ail 
temperatures  in  agreement  with  previous  flowing  afterglow  studies.'’^  The  n=1  cluster 
undergoes  ligand  switching  at  around  75%  of  the  collision  rate  at  all  temperatures.  The 
n=2  cluster  has  a  negative  temperature  dependence  in  the  rate  constant,  probably 
because  of  the  small  association  channel  present  at  150  K.  The  chemistry  of  the  n=3 
and  higher  clusters  is  more  complex.  The  rates  continue  decreasing  with  increasing 
cluster  size.  As  seen  in  Fig.  1  illustrates  where  the  reaction  efficiency  given  by  the  ratio 
of  the  observed  rate  constant  to  the  collision  rate  constant  is  plotted  as  a  function  of 
temperature  and  cluster  size,  the  largest  cluster  observed  at  each  temperature  shows 
an  anomalous  increase.  The  deviation  occurs  where  the  largest  cluster  is  only 
marginally  stable'll-  and  may  be  thermally  decomposing.  However,  this  would  affect 
the  rate  constant  of  the  next  smallest  cluster  which  is  unaffected  by  the  presence  of 
larger  clusters  in  the  flow  tube.  The  presence  of  the  association  channel  does  not 
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account  for  these  results  because  it  is  present  for  n=2  which  has  a  more  typical 
temperature  dependence  for  a  system  involving  an  endothermic  ligand  switching 
reaction  with  an  association  channel  at  low  temperature.  The  results  for  CH2O  show  that 
the  rate  constants  for  this  reaction  used  by  Yu  et  al.''^  are  too  large.l^ 

To  investigate  the  statistical  nature  of  Sn2  substitution  reactions  for  weakly 
ionized  plasmas,  the  reactions  Cl*  with  CaHsBr  and  n-CsHrBr  have  been  studied  as  a 
function  of  kinetic  energy  and  temperature  in  the  VT-SIFDT  from  220-500  K.  For  Cl 
reacting  with  C2H5Br.  the  Sn2  reaction  path  to  give  Br  dominates  at  all  temperatures 
and  kinetic  energies.  However,  an  association  channel  accounts  for  5%  of  the  products 
at  298  K  independent  of  kinetic  energy  and  disappears  at  500  K,  but  it  makes  a  larger 
contribution  at  220  K  showing  a  negative  kinetic  energy  dependence.  The  rate 
constants  also  show  a  negative  kinetic  energy  dependence  and  are  slightly  faster  than 
the  corresponding  pure  temperature  data.  Fig.  2  shows  the  total  and  Sn2  only  rate 
constants  as  a  function  of  kinetic  energy,  along  with  previous  experimental 
.  measurements. 20*22  The  rate  constant  for  the  Sn2  reaction  only  Is  obtained  by 

multiplying  the  Sn2  branching  fraction  by  the  total  rate  constant.  Comparing  the  data  at 
different  temperatures  with  the  same  kinetic  energy  reveals  that  the  internal 
temperature  does  not  have  a  strong  influence  on  the  rate  constants. 

Similar  trends  arise  for  the  reaction  of  Cl"  with  n-CsHjBr,  where  the  Sn2  channel 
again  predominates.  The  results  are  given  in  Fig.  3  along  with  a  previous  result.20  jh© 
rate  constants  are  essentially  independent  of  temperature  and  kinetic  energy  from  298 
to  500  K,  again  showing  a  5%  association  product  channel  only  at  298  K.  The  trend 
indicates  that  the  reaction  has  no  substantial  dependence  on  internal  energy.  However, 
a  factor  of  4  increase  in  the  rate  constant  at  zero  drift  field  at  230  K  occurs  and  the  rate 
constants  at  kinetic  energies  up  to  0.07  eV  decrease  with  increasing  kinetic  energy. 

The  higher  kinetic  energy  data  are  only  slightly  larger  than  the  values  at  higher 
temperatures.  The  association  channel  accounts  for  35%  of  the  products  at  zero  drift 
field  and  decreases  to  30%  at  0.07  eV  at  200  K.  A  second  mechanism  involving 
association  followed  by  dissociation  to  Bf  and  Cl"  via  collisions  with  He  buffer  gas  has 
been  examined  by  selectively  scavenging  the  association  product  via  secondary 
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reaction  with  HCI.  The  results  confirm  that  this  mechanism  occurs  in  addition  to  the 
direct  Sn2  reaction  for  n-CsHyBr. 

A  similar  test  with  CaHsBr  displays  analogous  results,  although  of  smaller 
magnitude.  The  energetics  of  the  reactions  with  both  alkylbromides  are  similar  as  well 
and  it  appears  that  the  importance  of  the  thermal  association  and  dissociation 
mechanism  depends  on  a  balance  between  the  well  depth,  barrier  height  and/or  the 
number  of  degrees  of  freedom.  Based  on  the  presence  of  the  association  channel, 
previous  work  with  Sn2  reactions  with  association  channels,23. 24  and  a  weak  internal 
and  translational  energy  dependence,  the  reactions  of  Cl'  with  CaHsBr  and  n-CsHyBr 
should  behave-statistically.25  This  finding  is  in  contrast  to  previous  work  with  CHaBr 
that  indicates  the  reaction  is  not  statistical. 26-29 

Two  reactions  that  have  been  studied  to  high  temperatures  are  with  CO2  and 
SO2,  which  are  minor  reactions  in  the  atmosphere.30-32  jhe  reaction  of  N2‘^  with  CO2 
has  been  studied  as  a  function  of  temperature  from  300-1400  K  in  the  HTFA  and  kinetic 
energy  from  0.04-0.27  eV  in  the  VT-SIFDT  at  298  K.  Fig.  4  shows  the  rate  constants 
plotted  vs.  the  average  total  energy  given  by  Eq.  2  for  a  temperature  T: 

{Total  Energy)  =  R£,„  +(K£)  +  (e*")  +(£"f"“>  (2) 

where  RE\on  is  the  recombination  energy  of  the  reactant  ion.  The  average  kinetic 
energy,  <K£),  is  determined  by  Eq.  1  in  the  VT-SIFDT  or  by  simply  3/2/cBrin  the  HTFA. 
The  average  internal  energies  of  the  reactant  ion,  {e^)  ,  and  neutral,  ,  include 

an  average  rotational  energy  of  n/2kBT  where  n  is  the  number  of  rotational  degrees  of 
freedom  and  an  ensemble  average  over  a  Boltzmann  distribution  of  vibrational  energy 
levels  for  all  modes  in  the  harmonic  approximation.  Zero  point  energy  is  neglected  in  all 
cases.  The  effective  temperature  in  the  drift  tube  experiments  used  to  evaluate  Eq.  2  is 
determined  by  Eq.  1  with  mn=mb. 

For  N2^  with  CO2,  the  rate  constant  at  298  K  with  no  drift  voltage  agrees  with 
those  from  the  HTFA  and  previous  experiments,33-36  proceeding  close  to  the  Langevin 
collision  rate  constant.  The  rate  constants  decrease  with  increasing  energy,  indicating 
the  reaction  proceeds  through  a  collision  complex  as  supported  by  previous  work.37, 38 
However,  the  HTFA  results  that  reflect  are  larger  fraction  of  internal  energy  are 
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consistently  lower  than  the  VT-SIFDT  results  with  most  translational  excitation.  Internal 
energy  thus  hinders  the  reactivity  relative  to  translational  energy.  These  results  contrast 
with  the  results  for  the  energetically  similar  Ar*  reaction  with  CO2  at  high  temperatures. 

In  that  experiment,  the  rate  constants  in  the  HTFA^s  and  a  flow  drift  tube  (FDT)‘^o  agree 
when  plotted  vs.  total  energy  up  to  0.4  eV.  Above  this  energy,  the  rate  constants  in  the 
HTFA  are  lower  than  the  FDT  values.  This  energy  corresponds  to  the  bond  energy  of 
the  (Ar*C02r  complex  which  may  interfere  with  the  charge  transfer  reaction.  The 
complex  should  have  similar  energy  with  N2'^:  however,  the  biggest  difference  is  the 
rotational  excitation  of  the  N2'^  ion  not  present  with  Ar^. 

The  reaction  of  N2^  with  SO2  has  also  been  examined  in  the  HTFA  up  to  1400  K. 
Fig.  5  shows  the  rate  constants  plotted  vs.  the  average  total  energy.  The  FDT  results  of 
Dotan  et  al.  are  also  shown,  converted  to  total  energy  using  an  effective  temperature 
derived  from  Eq.  1  as  described  above.  The  rate  constants  agree  within  the  error  of  the 
two  experiments  up  to  0.4  eV.  The  FDT  rate  constants  increase  above  0.6  eV  where  an 
endothermic  dissociative  charge  transfer  channel  to  give  SO^  is  accessible.  Fig.  6 
shows  the  rate  constants  for  the  individual  product  channels  in  both  experiments.  The 
rate  constant  for  generating  SO*  increases  rapidly  when  the  SO*  product  can  be 
formed,  but  the  trends  in  the  two  instruments  seem  to  differ.  However,  if  the  data  is 
replotted  as  the  ratio  of  the  SO*  channel’s  rate  constant  to  the  maximum  possible  rate 
constant  for  reaction  from  N2*  (v>0),  then  the  two  data  sets  concur  within  the  scatter. 
Consequently,  the  enhancement  in  the  rates  is  due  to  producing  SO*  from  vibrationally 
excited  states  of  N2*  as  seen  in  previous  work.^"*  With  Ar*  and  SO2,  the  vibrational 
excitation  in  SO2  at  elevated  temperatures  enhances  the  rate  constants;^^  however,  this 
enhancement  may  occur  with  N2*  in  a  region  where  the  effect  of  the  N2*(v^0)  states 

dominates.'^2 

To  address  the  needs  of  the  aforementioned  combustion  models,  the  kinetics  of 
the  reactions  of  air  plasma  ions  with  toluene,  ethylbenzene  and  propylbenzene  have 
been  studied  in  the  VT-SIFT  from  298-500  K.  Other  ions  with  recombination  energies 
between  9  and  16  eV  have  also  been  studied  at  298  K  in  the  VT-SIFT  to  construct 
breakdown  curves  for  the  dependence  of  the  branching  fractions  on  predominately 
electronic  energy.  These  curves  allow  comparisons  with  photodissociation 
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experiments.'^^  in  addition,  the  reactions  of  NO^  and  O2*  with  all  three  alkylbenzenes 
have  been  studied  up  to  1400  K  in  the  HTFA.  Comparing  the  results  from  the  two 
instruments  reveals  how  the  larger  fraction  of  internal  energy  in  the  high  temperature 
experiments  affects  the  reactivity  relative  to  the  mostly  electronic  excitation  in  the  VT- 
SIFT.  Table  3  shows  the  rate  constants  for  the  air  plasma  cation  reactions  at  298  K.  All 
reactions  proceed  at  the  Su-Chesnavich  collision  rate  constant  for  all  ions,44, 45 
including  those  not  shown  in  the  table.  The  rate  constants  for  the  NO"^  and  O2*  ions  are 
independent  of  temperature  to  1400  K,  consistent  with  reactions  that  proceed  at  the 
collision  rate.^- 

The  breakdown  curves  for  toluene,  ethylbenzene  and  propylbenzene  are  shown 
in  Figs.  8-10.  The  branching  ratios  are  plotted  against  the  average  total  energy 
calculated  using  Eq.  2.  To  examine  the  effect  of  the  charge-transfer  product  ion’s 
lifetime  on  the  branching  ratios,  the  flow  tube  studies  have  been  compared  to  guided- 
ion  beam  (GIB)  results  that  have  been  conducted  at  a  few  orders  of  magnitude  lower 
pressure.  Previous  experiments  have  shown  how  collisional  stabilization  of  a  long-lived 
energized  charge-transfer  complex  by  the  buffer  gas  in  the  flow  tube  can  affect  the 
observed  dissociation  thresholds, 9. 10  warranting  the  low  pressure  studies.  The  GIB  is 
equipped  with  a  high  temperature  octopole^®>  (HT8P)  that  allows  cross  sections  and 
branching  fractions  to  be  measured  as  a  function  of  both  temperature  (i.e.,  internal 
energy  content)  and  kinetic  energy  up  to  10  eV.  The  GIB  data  are  also  shown  in  Figs.  8- 
10  labeled  as  HT8P,  taken  at  298  and  650  K. 

Comparing  the  flow  tube  results  in  Fig.  8  for  toluene,  the  HTFA  and  SIFT  results 
agree  except  for  the  reaction  of  NO+  in  the  HTFA,  indicating  that  rotational  and 
vibrational  excitation  of  the  neutral  reagent  affects  the  reactivity  the  same  as  electronic 
energy.  However,  the  HT8P  results  indicate  that  translational  energy  is  not  as  efficient 
at  promoting  fragmentation.  The  NO^  data  show  a  fall-off  in  the  C/Hs^  branching 
fraction  at  the  thermochemical  threshold  for  C/Hr'^  formation.  Collisions  with  the  He 
buffer  gas  cause  dissociation  of  the  charge-transfer  product,  as  evidenced  by  a  lack  of 
such  behavior  at  similar  energies  in  the  HT8P  at  low  pressures.  However,  this  behavior 
has  been  observed  previously  in  the  flow  tube  with  CsHi/^  formed  from  isobutane.'^® 
Over  85%  of  the  dissociation  products  are  CtHt"^  in  equal  amounts  of  the  benzylium  and 
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tropylium  isomer  as  determined  in  a  previous  VT-SIFT  study.'^®  The  branching  ratios  for 
this  product  are  given  in  the  lower  panel  of  Fig.  8. 

The  observed  onset  of  dissociative  charge-transfer  with  toluene  in  the  flow  tube 
experiments  occurs  1  eV  higher  than  the  thermochemical  threshold.^o  jhe  CrHe'^  ion 
lifetime  is  such  that  an  extra  0.5  eV  of  energy  must  be  added  to  get  the  molecule  to  fall 
apart  on  the  timescale  for  observation  in  the  flow  tube,  creating  a  kinetic  shift.  However, 
the  threshold  is  shifted  another  0.5  eV  higher  because  collisions  with  the  He  buffer 
stabilize  the  energized  charge-transfer  product  ion.  Similar  trends  have  been  seen  in 
flow  tube  studies  of  benzene  and  naphthalene  to  1400  K.9- 

The  analogous  breakdown  curves  for  ethylbenzene  are  shown  in  Fig.  9.  The 
NO"^  reaction  in  the  HTFA  again  shows  the  unusual  trend  in  the  CeHio"^  branching 
fraction  caused  by  dissociation  induced  through  collisions  with  the  He  buffer.  The  fall- 
off  occurs  at  the  thermochemical  threshold  and  is  not  observed  in  the  HT8P  work  at 
similar  energy  but  lower  pressure.  The  O2  data  from  the  HTFA  agree  well  with  the  VT- 
SIFDT  data  intimating  that  rotational,  vibrational  and  electronic  energy  are  equivalent. 
Based  on  the  HT8P  results,  kinetic  energy  is  not  as  efficient  with  ethylbenzene. 
However,  it  is  more  efficient  than  it  is  for  toluene.  The  lower  panel  in  Fig.  9  shows  that 
over  90%  of  the  dissociation  products  are  C7H7‘^  of  which  70%  is  in  the  benzylium 
structure.49  The  better  agreement  between  the  low  kinetic  energy  data  from  the  HT8P 
and  the  flow  tube  data  at  1  Torr  pressure  indicates  that  less  of  a  kinetic  shift  and 
collisional  stabilization  occur.  This  trend  is  consistent  with  deuterium  isotope 
experiments  that  show  that  the  lifetime  of  the  energized  CaHio^  charge-transfer  product 
ion  is  shorter  than  that  for  toluene.5i 

Fig.  10  shows  the  breakdown  curves  for  propylbenzene.  As  with  the  other 
reactants,  the  charge-transfer  product  ions  from  reaction  in  the  HTFA  undergo 
dissociation  via  collisions  with  the  He  buffer  at  higher  temperatures.  Rovibrational  and 
electronic  energy  affect  the  reactivity  identically  and  kinetic  energy  is  the  most  efficient 
with  propylbenzene.  Approximately  90%  of  the  products  in  the  lower  panel  of  Fig.  10 
are  C7H7'^,  of  which  90%  is  in  the  benzylium  form.49  The  flow  tube  data  at  1  Torr 
pressure  and  HT8P  kinetic  energy  data  at  10“^  Torr  have  even  better  agreement  than 
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with  ethylbenzene  at  low  energies,  consistent  with  the  shorter  lifetime  of  the  C9H12'' 
energized  ion.^^ 

The  data  from  the  combined  experiments  provides  insight  into  the  charge- 
transfer  dynamics  for  alkylbenzene  reactions  with  air  plasma  ions.  At  low  kinetic 
energies,  the  rate  constants  proceed  at  the  collision  rate.  Also,  internal  and  electronic 
energy  behave  similarly,  together  indicating  that  complex  formation  occurs.  However, 
kinetic  energy  is  not  equivalent  to  the  other  forms  of  energy:  therefore,  complex 
formation  does  not  completely  describe  the  mechanism.  The  reduced  efficiency  of 
kinetic  energy  implies  that  long-range  charge-transfer  processes  occur  at  distances 
whereJittle  kinetic  energy  is  exchanged.  At  higher  kinetic  energy,  the  threshold 
behavior  of  NO"^  reacting  with  ethylbenzene  and  propylbenzene  in  the  HT8P  is  similar  to 
line-of-centers  functions  that  describe  collision  induced  dissociation  (CID)  cross 
sections.53  These  CID  cross  sections  peak  at  a  value  approaching  the  effective  hard 
sphere  cross  section.  The  different  maximum  amounts  of  fragmentation  with  and 
NO^  in  the  HT8P  may  arise  from  differences  in  the  hard  sphere  cross  sections. 
However,  the  asymptotic  amounts  of  dissociation  also  depend  on  the  initial  total  energy 
available.  Therefore,  near-resonant  processes  may  also  occur.  1 1 

The  ion-molecule  chemistry  of  weakly  charged  plasmas  and  high  temperature 
environments  has  been  explored  using  the  fast  flow  tube  methods  available  at  the 
Space  Vehicles  Directorate  of  AFRL  at  Hanscom  AFB.  The  unique  capabilities  provide 
comprehensive  coverage  over  the  range  of  operating  temperatures  found  in  Air  Force 
systems.  A  temperature  variable  version  of  a  recently  constructed  turbulent  ion  flow 
tube54  (TIFT)  will  expand  these  capabilities  to  afford  examining  the  effects  of  pressure 
and  temperature  on  ion-molecule  chemistry.  It  is  this  adaptability  which  will  allow  these 
powerful  techniques  to  address  critical  issues  for  the  understanding  of  basic  processes 
influencing  crucial  Air  Force  systems. 
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Table  1.  Temperature  dependent  rate  constants  and  collision  rate  constants  for  the 
reactions  of  H30'^(H20)n  with  acetaldehyde. 

Rate  Constant  (10‘®  cm^  s'^),  [kd 

Reaction  200  K  294  K 

HaO^  +  CH3CHO  CHaOH^  +  H2O 

H30^(H20)  +  CH3CHO  ^  H30^(CH20)  +  H2O 

H30^(H20)2  +  CH3CHO  H30*(H20)(CH3CH0)  +  H2O 
H30^(H20)3  +  CH3CHO  ->  Products 
H30^(H20)4  +  CH3CHO  Products 


4.7  [4.7]  3.7  [3.9] 

4.1  [3.8]  3.5  [3.2] 

3.7  [3.5]  3.1  [2.9] 

3.4  [3.3]  2.9  [2.7] 

3.5  [3.1] 
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Table  2.  Temperature  dependent  rate  constants,  collision  rate  constants,  and  product 
distributions  for  the  reactions  of  H30'^(H20)n  with  formaldehyde. 


Rate  Constant  (10  ®  cm®  s*^),  [kd 
Branching  Fraction 


Reaction 

150  K 

200  K 

294  K 

HaO^  +  CH2O  ^ 

Products 

4.4  [4.9] 

3.8  [4.3] 

3.3  [3.6] 

CH20H^  +  H2O 

1.0 

1.0 

1.0 

H30^(H20)  +  CH2O  -> 

Products 

3.1  [4.1] 

2.8  [3.6] 

2.3  [3.0] 

H30^(CH20)  +  H2O 

1.0 

1.0 

1.0 

H30*(H20)2  +  CH2O 

Products 

2.0  [3.8] 

1.6  [3.3] 

1.1  [2.8] 

H30*(H20)(CH20)  +  H2O 

-0.5 

1.0 

1.0 

H30^(H20)2(CH20) 

-0.5 

H30^(H20)3  +  CH2O 

Products 

1.6  [3.6] 

1.2  [3.2] 

1.4  [2.7] 

H30^(H20)2(CH20)  +  H2O 

-0.5 

1.0 

1.0 

H30^(H20)3(CH20) 

-0.5 

H30^(H20)4  +  CH2O 

Products 

0.7  [3.5] 

1.7[3.1] 

H30*(H20)5  +  CH2O 

Products 

1.1  [3.5] 
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Table  3.  Rate  constants  in  cm3  s-1  for  the  reaction  of  air  plasma  ions  with 
alkylbenzenes  at  298  K  as  measured  in  the  VT-SIFDT.  The  collision  rate  constants 
have  been  calculated  using  the  Su-Chesnavich  model  ^5 


Rate  Constant,  [Collision  Rate]  (10  ®  cm®  s*^) 

Ion 

Toluene 

Ethylbenzene 

Propylbenzene 

(Recombination 

Energy) 

CrHs 

CsHio 

C9H12 

NO^ 

(9.26  eV) 

1.6  [1.8] 

2.0  [2.0] 

2.1  [2.1] 

O2" 

(12.07  eV) 

1.7  [1.8] 

2.2  [1.9] 

1.9  [2.0] 

0^ 

(13.62  eV) 

2.2  [2.3] 

2.4  [2.6] 

2.6  [2.7] 

(14.53  eV) 

2.2  [2.5] 

2.6  [2.7] 

2.7  [2.8] 

N2^ 

(15.58  eV) 

1.9  [1.9] 

1.9  [2.0] 

2.1  [2.1] 
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Fig.  1.  (a)  Reaction 
H30^(H20)n  with  Ch 
200  K,  and  294  K 


Fig.  2.  Rate  constants  in  cm^  s"^  for  the  reaction  of  Cl’  with  CaHsBr  as  a  function  of 
kinetic  energy  in  eV.  Closed  and  open  circles  represent  the  total  rate  constant  and  the 
Sn2  only  rate  constant  at  220  K,  respectively.  Closed  and  open  squares  represent  the 
total  rate  constant  and  the  Sn2  only  rate  constant  at  298  K,  respectively.  The  pure 
temperature  dependent  data  of  the  direct  Sn2  reaction  is  given  by  triangles  with  a  least 
squares  fit  line.  The  triangles  include  data  at  400  and  500  K  where  no  correction  is 
needed.  For  the  correction  at  low  temperature,  see  the  text.  The  solid  line  refers  to  the 
study  of  Caldwell  et.  al.  (CMK).21  The  dashed  line  refers  to  the  study  of  Knighton  et  al. 
(KBOG).22  jhe  X  refers  to  the  study  of  Depuy  et  al.  (DGMB).20 


Cr  +  CH  Br->Br +CH  Cl 

2  5  2  5 


Cr(C  H  Br) 


0.02  0.03  0.04  0.05  0.06  0.07  0.08 


Kinetic  Energy  (eV) 
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Fig.  3.  Rate  constants  in  cm^  s'^  for  the  reaction  of  Cl’  with  n-CsHrBr  as  a  function  of 
kinetic  energy  in  eV.  Closed  and  open  circles  represent  the  total  rate  constant  and  the 
Sn2  only  rate  constant  at  230  K,  respectively.  Closed  and  open  squares  represent  the 
total  rate  constant  and  the  Sn2  only  rate  constant  at  298  K,  respectively.  The  diamonds 
refer  to  500  K  data.  The  x  refers  to  the  study  of  Depuy  et  al.  (DGMB).20  The  point 
labeled  “Direct  only”  is  an  upper  limit  to  the  direct  Sn2  reaction  at  230  K.  See  text  for 
details. 


Cr  +  CH  Br->Br +CH  Cl 

3  7  3  7 


0.02  0.03  0.04  0.05  0.06  0.07  0,08 

Kinetic  Energy  (eV) 
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Fig.  4.  Rate  constants  in  cm^  s'^  for  the  reaction  of  N2''  with  CO2  plotted  as  a  function  of 
the  average  total  energy  in  eV.  The  solid  symbols  (► )  represent  the  current  HTFA 
results  and  the  open  symbols  (□)  represent  the  current  SIFDT  results  at  298  K.  The 
dashed  line  represents  the  Langevin  collision  rate  constant. 
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Fig.  5.  Rate  constants  in  cm^  s'^  for  the  reaction  of  N2^  with  SO2  plotted  as  a  function  of 
the  average  total  energy  in  eV.  The  solid  symbols  (►  )  represent  the  current  HTFA 
results  and  the  open  symbols  (□)  represent  the  previous  FDT  results  of  Dotan  et  al. 

The  dashed  line  represents  the  Su-Chesnavich  collision  rate  constant.44. 45  jhe  HTFA 
rate  constants  have  been  multiplied  by  1 .2  to  better  overlap  with  the  FDT  data.  The 
solid  lines  are  fits  to  the  data  using  the  functional  form:  k(E)  =  Ai(0.1 1 1/E)"  +  A2  exp(- 
Ea/E). 
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Fig.  6.  Rate  constants  in  cm^  s'^  for  30“^  production  only  for  the  reaction  of  N2^  +  SO2 
plotted  as  a  function  of  the  average  total  energy  in  eV.  The  solid  symbols  (►  )  and  (■) 
represent  the  current  HTFA  results  for  the  SO2*  and  SO"^  channels,  respectively  and  the 
open  symbols  (O)  and  (□)  represent  the  previous  FDT  results  of  Dotan  et  al.55  for  the 
302'^  and  30'"  channels,  respectively. 
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Fig.  7.  Rate  constants  for  the  SO^  channel  divided  by  the  maximum  rate  constant  for 
the  reaction  of  N2'^  (v  ^ ).  The  (v  ^ )  rate  is  defined  as  the  collision  rate  constant 

times  the  N2^  (v^)  population.  The  solid  symbols  (► )  represent  the  current  HTFA 
results  and  the  open  symbols  (□)  represent  the  previous  FDT  results  of  Dotan  et  al.55 
The  open  triangles  (A)  represent  the  FDT  with  a  background  subtracted. 
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Fig.  8  Breakdown  diagrams  for  the  reactions  of  toluene  (CrHs)  with  various  ions 
showing  (a)  and  (b)  CrHr'^  as  a  function  of  total  energy.  An  <8)  symbol  marks  an 
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O2*  SIFDT  result  at  300  K  with  a  center-of-mass  collision  energy  of  0.52  eV 
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Fig.  9  Breakdown  diagrams  for  the  reactions  of  ethylbenzene  (CsHio)  with  various  ions 
showing  (a)  CsHio"^  and  (b)  as  a  function  of  total  energy.  An  <S)  symbol  marks  an 
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SIFDT  result  at  300  K  with  a  center-of-mass  collision  energy  of  0.52  eV 
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Fig.  10  Breakdown  diagrams  for  the  reactions  of  propylbenzene  (C9H12)  with  various 
ions  showing  (a)  C9Hi2^  and  (b)  as  a  function  of  total  energy.  An  <S>  symbol  marks 
an  02"^  SIFDT  result  at  300  K  with  a  center-of-mass  collision  energy  of  0.52  eV 
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